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ABSTRACT 


Experimental work in support of stress studies in high speed 
silicon sheet growth has been emphasized in this quarter. Creep 
experiments utilizing four-point bending have been made in the 
temperature range from 1000°C to 1360°C in CZ silicon as well as on EFG 
ribbon. Creep is found to be strong, even at the lower temperatures. 
Important differences in primary creep are observed with silicon 
impurity content in the higher temperature regions. 

A method to measure residual stress over large areas using laser 
interferometry to map strain distributions under load is under 
development. The experimental system to make these measurements on CZ 
silicon has been built and tested. 

A fiber optics sensor to measure ribbon temperature profiles has 
been constructed and is being tested in a ribbon growth furnace 
environment. It has been necessary to enclose the optical fibers in a 
water-cooled jacket to improve reproducibility and reliability of the 
sensor . 

Stress and temperature field modeling work has been directed 
toward improving various aspects of the finite element computing 
schemes. Difficulties in computing stress distributions with a very 
high creep intensity and with non-zero interface stress have been 
encountered and additional development of the numerical schemes to cope 
with these problems is required. Temperature field modeling has been 
extended to include the study of heat transfer effects in the die and 
meniscus regions and efforts are underway to improve the overall scheme 
of calculating interface region temperature profiles for real system 
geometries. 


"This report was prepared as an account of work sponsored by the United 
States Government. Neither the United States nor the United States 
Department of Energy, nor any of their employees, nor any of their 
contractors, subcontractors, or their employees, makes any warranty 
express or implied, or assumes any legal liability or responsibility 
for the accuracy, completeness or usefulness of any information, 
apparatus, product or process disclosed, or represents that its use 
would not infringe privately owned rights." 
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I . INTRODUCTION 


A satisfactory model that can account for stresses generated in silicon 
sheet grown at high speeds is not yet available. Many attempts to calculate 
residual stresses have been made, but all of these suffer from inadequacies in 
one area or another. This report describes the work in progress under this 
subcontract to attempt to develop and test in the laboratory a stress-temperature 
field model for silicon ribbon EFG. In one subtask, a computer code developed at 
Harvard University to calculate stress with plastic deformation in steady-state 
sheet growth is being applied to study temperature field-stress relationships. 
The stress state is parameterized by a two-dimensional temperature field and 
growth speed. Incorporation of time dependent stress relaxation effects is 
through a creep law to model the impact of plastic flow on the sheet residual 
stress state. A second aspect of the program deals with the development of a 
model to predict the temperature field in a moving sheet from given system 
component temperatures (i.e., the sheet environment), and studies experimental 
means to verify the model. This work is attempting to integrate temperature 
field modeling in the solid which takes into account real system geometries with 
a more detailed heat transfer model developed at MIT for the die top and meniscus 
regions below the growth interface. 

The predictions of the stress analysis and temperature field models are 
supported on a qualitative level by experimental results obtained from 
characterization of 10 cm wide ribbon grown in an EFG cartridge system. Residual 
stresses in both buckled and unbuckled as-grown ribbon are much lower than 
predicted for the thermoelastic case. These can only be achieved through high 
rates of creep for the temperature profiles that are calculated for the 10 cm 
ribbon. Defects are observed in the ribbon cross section which originate from 
stress instabilities. 

Both modeling and experimental work are in progress to attempt to apply 
the stress analysis at a more quantitative level. The modeling is examining 
stress-temperature field relationships in a search for low stress growth 
configurations. The experimental efforts are in study of the high temperature 
creep response of silicon to provide information on the constitutive 
relationship, development of a technique to measure the temperature field in the 
sheet during growth, and of means to measure residual stress distributions. A 
temperature sensor based on fiber optics has been constructed for the temperature 
field measurements. A laser interferometry technique is being evaluated for 
residual stress measurements at the University of Illinois. 

Attempts to verify aspects of the stress and temperature field modeling 
with the 10 cm cartridge system have met with limited success. Due to the 
complexity of the growth interface environment it has not been possible to use 
the fiber optics sensor. The cartridge system has lacked flexibility to allow 
imposition of changes in growth conditions and system temperature fields in a 
known and controllable manner in order to study temperature field-stress 
relationships at a quantitative level. A system that has the required 
flexibility for testing of the stress analysis model is being constructed. 
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II. PROGRESS REPORT 


The work carried out on development of understanding of stress generating 
mechanisms in vertical silicon sheet growth is described in the next four 
sections . 


A. Advanced System Design 

1. Temperature Field Modeling (H.M. Ettouney and R.A. Brown, MIT) 

The finite element analysis techniques developed for solution of 
coupled capillary and heat and mass transfer equations in the EFG system have 
been applied to study the 10 cm cartridge system [1]. Experimental data have 
been obtained for known environment temperature configurations of the ribbon on 
the thickness-velocity relationship and on aluminum segregation through the 
ribbon thickness. The latter gives information on the interface shape. The data 
will be used to test certain aspects of the heat transfer model in use and thus 
to obtain information on the dominant heat transfer parameters influencing ribbon 
temperature distributions. This information will be used in stress analysis to 
help in placing these calculations on a more quantitative level. 

In the following, we report the fitting results of the theoretical 
predictions by the finite element model developed in [1] and the experimentally 
measured field variables. In the first part, variations in the crystal thickness 
as a function of growth rate are fitted by solution of the thermal /capil lary 
problem, which includes simultaneous solution of the energy balance equations in 
melt, solid, and die and the Young-Laplace stress balance equation. In the 
second part, dopant segregation along the melt/solid interface is compared to the 
experimental measurements. The theoretical results are obtained by solving the 
momentum and the mass transfer equations using the technique presented in [2]. 


a . Fitting of the Crystal Thickness - Velocity Data 

Figure la shows the measured variations in the crystal 
thickness as the growth rate is increased from 2.5 to 4 cm/minute. Figure lb 
shows the axial temperature profile of the surroundings during the experiment. 
Crystal thickness measured at the right, center and left sides of the sheets is 
shown as curves a, b, and c in Fig. la, respectively. The left and right 
measurements are made at a distance of 1 cm from the edges of the crystal. The 
changes in the crystal thickness along the same horizontal level are caused by 
nonuniformities in the rate of heating in the die. The surrounding temperature 
profile shown in Figure lb is driven by a water cooler and afterheater. The 
cooler is placed at a distance of 0.2 cm from the top of the die and has a length 
of 0.3 cm. The afterheater has a constant temperature gradient of 50 K/cm and 
its high temperature setting is equal to 1473 U K. At both ends of the cooler 
there are rapid increases in the surrounding temperature followed by a linear 
increase to the heater temperature, downstream of the cooler, and by a linear 
increase to a temperature close to the die temperature upstream of the cooler. 

The dimensionless temperatures of the surroundings are shown on Figure lb. 
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The finite element model described in detail in reference [1] is used to 
fit curve a in Figure la. In this model the values of the thermophysical 
properties used to characterize molten and solid silicon and the die material are 
the same as those given in Table II of reference [1], except for the thermal 
conductivity of solid silicon which is taken to be temperature independent there. 
Here, it is redefined according to the relation K(x,y) = (K^/T„)/T(x,y^ , where K M 
is the solid phase thermal conductivity at the melting point, 0.22 J/s°K/cm. The 
height and length of the die base and tip are the same as those given in 
reference [1], Table I. The capillary spacing and the static head are taken 
equal to be 0.05 cm and 4 cm, respectively, and the surrounding temperature 
profile given in Figure lb is used. These specifications define all the system 
parameters except for the die base temperature, T , and the surrounding 
temperature at the die base, T^,. The two unspecified parameters are carefully 
varied within a realistic range; i.e., T^ + 20 < T < T M + 50°K and T° ^ T , in 
order to obtain the best theoretical fit to the experimental measurements. The 
resulting fit is shown in Figure 2 and the calculated axial temperature profile 
along the surface at a growth rate of 4 cm/minute is shown together with the 
surrounding temperature profile in Figure lb. The best fit was obtained for T = 
1727°K and T° = 1733°K. The sensitivity of these results with respect to 0 
variations in the two fitting parameters was examined. It is found that to 
achieve a 50% increase in the crystal thickness, the die temperature needed to be 
decreased by 13°K (compare the two curves in Figure 2 for T = 1727°K and T = 
1714°K). On the other hand, large adjustments in T^, at a constant T* , were 
required to achieve the same degree of change in the crystal thickness. For 
example, the crystal thickness was increased from 0.008 to 0.013 cm, at a growth 
rate of 4 cm/minute and a die temperature of 1727°K, when T^, was decreased from 
1961°K to 1372°K . 


b . Fitting of Dopant Segregation 

Calculations of the dopant distribution in the melt and its 
segregation along the melt/solid interface are more complicated than the heat 
transfer calculations since they require calculations of the flow field. The 
mass and momentum equations and the finite element approximation for the field 
variables are given in detail in [2]. The fit for the dopant segregation is 
compared with the experimentally measured aluminum profiles shown in Figure 3. 

The three profiles are measured along the crystal thickness at right, center and 
left, Figures 3a, 3b and 3d, of the crystal at a growth rate of 2.6 cm/minute for 
the same samples as in Figure la. In these three profiles the average aluminum 
segregation and operating conditions of the solid curve in Figure 2, defined as 
the ratio between the difference of the solute concentration at the surface and 
the centerline to the bulk concentration, varies between 10 and 20%, excluding 
that caused by a structural defect (e.g., a grain boundary) in the center of 
Figure 3c. 

To analyze the theoretical predictions of the solute profiles, one has to 
take into account shapes of the melt/solid interface and the flow field pattern. 
In Figure 4a, the displayed melt/solid interfaces are concave, as viewed from the 
melt, and their deflections, defined as the ratio between the difference of the 
interface height at the centerline and the surface to the crystal thickness, are 
1.1, 1, 0.9. and 0.78 at growth rates of 2.5, 3, 3.5, and 4 cm/minute. As 
discussed in references [2] and [3], such interface shape always results in a 
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Comparison of experimental and calculated t-V 
relationships for 10 cm wide ribbon. 



crc* 



Fig. 3. Spreading resistance (resistivity) profiles for aluminum-doped 

ribbon grown at 2.5 cm/min. Traces (a), (b), and (c) are for the 
ribbon represented in curves (a), (b), (c) in Fig. la. The traces 
were taken down a 2°52' level through the ribbon thickness. 
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Interface shapes (a) and associated die top and meniscus melt streamlines (b)-(e) 
calculated for ribbon growth in system represented in Fig. 2. 




solute concentration that is higher at the centerline than at the crystal 
surface. Figures 4b, 4c, 4d, and 4e show the calculated streamlines at growth 
rates of 2.5, 3, 3.5, and 4 cm/minute. In all cases, a converging flow field 
develops. As a result, higher solute concentration is also obtained at the 
centerline, because of lateral convection. 


The isoconcentration lines and the interface concentration displayed in 
Figure 5 are calculated for a solute with equilibrium distribution coefficient of 
0.1 and a diffusion coefficient equal to that of aluminum (0.0005 cm /s) [4]. As 
shown in Figures 5b, 5c, 5d, and 5e, the isoconcentration lines for the 
interfaces in Figure 5a become more squeezed into the melt/solid interface as the 
growth rate is increased. Regardless of the increase in the growth rate, 
magnitude of the dopant segregation along the interface is decreased from 10.56, 
9.27, 7.84 to 6.64 at growth rates of 2.5, 3, 3.5, and 4 cm/minute, respectively. 
This behavior is caused by the continuous decrease in the net deflection of the 
melt/solid interface as the growth rate is increased. Further decreases in the 
magnitude of the equilibrium distribution coefficient results in a very small 
change in the segregation magnitude; e.g., at a growth rate of 2.5 cm/minute 
dopant segregation is equal to 10.56, 11.29, and 11.54 for K of 0.1, 0.05, and 
0.01, respectively. 


These results suggest that for aluminum, K = 0.002; segregation as 
predicted by the finite element model will remain within the range measured 
experimentally. 


The success of the initial fitting between the predictions of the finite 
element model and the experimental measurements suggest further application of 
the model is useful to simulate and study other EFG configurations and also to 
continue analysis of the available experimental data. In the next step, the 
Marangoni convection will be invoked in both of the thermal and mass transfer 
models. The results of this step should confirm the absence of Marangoni 
convection in the experiment, since very large solute segregations are expected 
to occur because of the high intensity of the flow field. These calculations 
will be followed by investigation of the system stability caused by 
thermo cap illary effects, where at high growth rates and static heads the meniscus 
becomes unstable and as a result the crystal pulls free. This behavior has been 
modeled and studied for a prototype EFG configuration (using the same system 
parameters of reference [1]), where Newton's method and the continuation 
technique were employed to calculate the critical points of stability. 

2. Stress Analysis (J.C. Lambropoulos and J.W. Hutchinson, Harvard 

University) 

The calculations of stress distributions with high rates of 
plastic deformation were reported in the last quarterly [5] to be only 
approximately correct. More detailed investigation into the numerical schemes 
that are used in the analysis has been carried out to check various aspects of 
the calculation. The cause of the difficulties in obtaining physically 
meaningful solutions has been traced to an inadequacy of the grid size very close 
to the growth interface to handle the very large changes in stress and strain 
rate that are imposed for the high creep intensity predicted by the data of 
Siethoff and Schroter [6], 
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for ribbon growth in system represented in Fig. 



The difficulty in dealing with the high rate of change of the variables 

near the growth interface appears to also occur in modeling of non-zero interface 

stresses. Rather than attempt to increase grid size to deal with this problem, 

investigation of a simpler scheme to study the high strain rate cases has been 

initiated. It appears that a physically valid approximation in this situation is 

to assume that e , the transverse shear strain rate, is essentially zero very 
VY • • J J 

close to the growth interface. This leads to some simple relations between the 

creep and thermoelastic strain rate components acting on the sheet. It is also 

evident that in this limit, the use of a zero interface stress component a 

boundary condition is justified. Essentially, the calculation results theX^ 

become independent of cr 

yy 

Details of this new approximation to account for high interface stresses 
will be reported in the future. 

B. Creep Law Studies 

Creep response of silicon has been studied in four-point bending of 
both CZ silicon and of ribbon at temperatures above 1000°C. The results on the 
defect structure are reported below. 

1 . CZ Silicon 

'< CZ silicon wafers with (100) surface orientations were subjected 
to incremental loading at 1215°C and at 1360°C. Th^ loads were chosen to arrive 
at comparable total strains of the order of 1 x 10 

The loading conditions were described in the previous quarterly report 
[5] . The stress acting on the samples was calculated assuming that a linear 
stress distribution exists across the thickness of the samples, as shown in Fig. 
6. Loads were applied by increasing the force on the sample with four-point 
bending until maximum sample surface stresses were of the order of 10 MPa. After 
a final displacement of the order of 0.1 cm was obtained, the load was removed 
and the samples cooled to room temperature. 

Cross-sectional micrographs of the CZ samples have been made along a 
longitudinal plane perpendicular to the bending axis, as illustrated in Fig. 6. 
The dislocation structures observed at locations of the samples marked A, B, and 
C in Fig. 6 are given in Figs. 7 and 8 for the 1215 C and 1360 C experiments, 
respectively . 

At 1215°C, the unstressed region A of the sample remains dislocation free 
(Fig. 7(a)). The high density of etch pits appears to be related to B-swirl 
defects found often in as-grown CZ silicon. This contr^sts^to the results at 
1360 C, where a dislocation density in excess of 1 x 10 /cm is produced. It is 
possible that these dislocations have migrated into the stress-free region, given 
the long time (3.5 hours) during which the sample was under stress. A second 
source of dislocations is from point defect condensation, whereby the B-swirls 
are assumed to coalesce into A-swirls of dislocation loop-type defects. 
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Region C, Np = 1.4E6/cm . 

(100) Czochralski silicon stressed at 1 21 5°C. 







In the transition region B of intermediate^str^sses , samples at both 
temperatures show dislocation densities of 1 x 10°/cm and above (Figs. 7(b) and 
8(b)). In addition, spatial nonuniformities occur, possibly due to the onset of 
polygonization . 

In the region C of highest, and most uniform stress, dislocation densities 
are at 1 .4 x 10°/cin and 1.8 x 10 /cm for the 1215°C and 1360°C samples, 
respectively. Again, distinct banding that leads to nonuniform dislocation 
distributions is observed. 

Some differences in the distributions through the thickness are observed 
at the two temperatures. This may be due to orientational effects as no effort 
was made to orient the samples with respect to the bending axis in these initial 
experiments. Future experiments with oriented wafers are planned. 

2 . Ribbon 


Ribbon was stressed in the four-point bending apparatus at 1050°C. 
The purpose of these experiments was to examine the extent of the creep response 
that would arise under relatively short times of loading. The loading conditions 
and displacements are given in Table I. After the sample cooled, the remaining 
deflection was measured and compared to the deflection under load. The results 
given in Table I show that the creep in the loading time of 40 minutes is 
sufficient to make the strain permanent in most cases, and relatively little 
elastic spring-back occurs. Considerable variability in this data was observed 
in sample to sample, as is seen in the representative displacement-time graphs 
shown in Figs. 9(a) and 9(b) for some samples. This is attributed to 
inhomogeneities in the material itself as well as to geometric nonuniformities 
which result in nonreproducib le loads. Applied stresses were in the range of 
30-35 MPa for the light loads, and 80-90 MPa for the heavier ones. 

It can be concluded from these experiments that creep can be still a very 
important factor in stress relaxation at temperatures of the order of 1000 C 
during ribbon growth. The above results are for as-grown ribbon samples which 
have already undergone considerable work hardening in the in situ furnace cool. 
Yet they still have a fast transient primary creep response under load. The data 
in Fig. 9 shows this is over in a time of the order of a minute or so. After 
steady-state creep sets in, it produces nearly complete relaxation in 20 minutes 
after application of the second load. 

The samples were further cross-sectioned and the defect structure 
examined. It has been difficult to establish the nature of the dislocation 
structures caused by the additional creep in view of the already high dislocation 
density in the as-grown material. Additional work is being done to look at these 
samples in more detail. 

3 . Transient Creep Response 

Emphasis here is on relatively short time (_< 10 minutes) 
non-steady-state behavior. Displacement time graphs have been generated for 
creep response of CZ and FZ silicon in order to compare the^ffect of impurities. 
The CZ wafer had measured substitutional carbon at 4.5 x 10 at/cc and oxygen at 
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Table I. Four-Point Bending at 1050 C. 
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Fig. 9. Response of ribbon stressed in four-point bending 
at 1050°C under loads of L, = 170g and L 9 = 440g. 
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1.0 x 10 at/cc. The FZ wafer had nondetectable levels ( <^ 1 x 10 at/cc) for 
these impur it ies . 

The displacement-time graphs generated at 1215°C for 10-second bending of 
the FZ (curve 1) and CZ (curve 2) and for 10 minutes for the CZ (curve 3) are 
shown in Fig. 10. It is evident that there is an order of magnitude difference 
in the strain rate between the FZ and CZ silicon, which can be attributed to 
impurities. The strain and strain rate data will be examined in relation to the 
defect distributions generated in these samples. 

C . Fiber Optics Temperature Sensor 

In the initial tests of the fiber optics based temperature sensor in 
the 10 cm cartridge, reliable and reproducible results were not obtained. Large 
variations in signal levels and differences from the calibration results obtained 
in a small test furnace occurred. The origin of the stray radiation was traced 
to the cartridge afterheater region, where the fibers passed very close to the 
active afterheater, rather than the probe tip. This was shown when the probe was 
rebuilt so that the tip was imbedded in the cold shoe and kept at temperatures 
below 300°C. Even then, the output signal was found to be dominated by 
variations in the afterheater temperature setting. 

After these initial attempts to install the optical fiber sensor in the 
cartridge failed, it was decided to construct a new optical fiber probe. The 
entire optical fiber section will be enclosed in a water-cooled jacket, which 
will extend to within several mm of the probe tip. This probe will not have the 
flexibility of the initial one, which consisted essentially of the fibers 
themselves. An additional difficulty at this time arises because only short 
lengths of thick (1 mm) Al^O^ fibers are available, rather than the longer 
thinner (0.25 mm) ones desired. This requires use of quartz fibers and completes 
the construction and assembly of the entire probe. Some problems with alignment 
at the terminations of the Al^O, and quartz have been experienced. We are still 
investigating the possibility of obtaining the longer A^O^ fibers. 

D . Residual Stress Measurements 


A residual stress measurement technique using laser interferometry to 
measure sample deflection under load is being evaluated at the University of 
Illinois. The technique has been successfully applied to calculate applied 
in-plane stresses in CZ wafers. A description of the technique and results of 
this preliminary study are contained in Appendix III. 

Modifications are under study to allow application of this method in the 
case of EFG ribbon. The main problem areas center around dealing with the 
irregular dimensions (surface striations and thickness nonuniformities) of the 
ribbon, which increases in complexity the fringe pattern generated, and in 
changes in fixturing and loading that are needed to accommodate the rectangular, 
rather than circular geometry. 
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APPENDIX II 


The following papers were presented at the Research Forum on High Speed 
Growth and Characterization of Crystals for Solar Cells, Port St. Lucie, Florida, 
July 25-27, 1983. 

a. "The Effect of Creep on the Residual Stresses Generated During Silicon 
Sheet Growth", by J.W. Hutchinson and J.C. Lambropoulos , Harvard 
University. 

b. "Mathematical Modeling of High-Speed Ribbon Systems: A Case Study of 

Edge-defined Film-fed Growth", by H.M. Ettouney and R.A. Brown, MIT. 

c. "Experimental Aspects of the Study of Stress Generating Mechanisms in 
Silicon Sheet Growth", by J.P. Kalejs and R.O. Bell, Mobil Solar' 
Energy Corporation. 

They will appear in the proceedings of the Forum to be published by JPL. 
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APPENDIX III 


Paper by A.T. Andonian and S. Danyluk on "Non-Destructive 
Determination of Residual Stresses in Circular Silicon 
Wafers", submitted for publication. 


Preceding Page Blank 


27 



NON-DESTRUCTIVE DETERMINATION OF RESIDUAL STRESSES 
IN CIRCULAR SILICON WAFERS 


A. T. Andonian and S. Danyluk 
Department of Civil Engineering, Mechanics and Metallurgy 
University of Illinois at Chicago 


Introduction 


In many instances, a predominant factor contributing to the failure of 
structural components is the magnitude of residual stresses. Generally, 
restraint limits the distortion but in some instances large local 
strains can arise which increase the risk of failure. Residual stresses 
can also contribute to fatigue failure, brittle fracture, structural 
instability, stress corrosion and in some cases cause premature buckling 
and warping or result in a dramatic reduction in the components 
resistance to brittle fracture. 


Manufacturing and fabricating processes are the most common causes of 
residual stress. Processes such as casting, welding, cutting, 
machining, molding and heat treatment introduce residual stresses [1-6]. 
For example, residual stresses can arise from thermal gradients in the 
components in processing and the unavoidable restraint offered by the 
surrounding structure. In addition, stresses can be generated during 
metal shrinkage. These stresses can be minimized to some extent by a 
careful control of the cooling rate, material chemistry or 
transformation characteristics, and both the internal and external 
restrain conditions. 


Residual stresses are particularly important in non-me tall ic plates such 
as semiconductor or ceramic sheet materials because the geometries are 
particularly conducive to buckling which in some cases can lead to 
fracture. For example, in the semiconductor industry, silicon is 
produced in thin sheet for the manufacture of large scale integrated 
circuits and solar arrays. In these cases growth conditions may 
influence precipitate (oxide) formation which leads to buckling of the 
ribbon in addition to influencing the electrical properties. 
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Measurement of Residual Stresses 

Measurement of residual stress cannot be accomplished by the 
conventional non-destructive methods of experimental stress analysis. 
The strain sensor (strain gage, photoelastic coating. Moire grating, 
etc.) being totally insensitive to the loading history of the component, 
measures only changes in strain after installation of the sensor. 

There are several non-destructive techniques available for determining 
residual stresses [7-13]. The two most widely used are the 'X-ray 
diffraction technique' and the 'blind-hole drilling' method. Both 
techniques have limited applicability depending on the material and 
component size. Aside from the bulk and complexity of the equipment, 
the X-ray technique is limited by the X-ray absorption which occurs 
within several micrometers of the surface depths and limits the 
measurement of residual stress to these depths. The resolution of this 
technique is approximately 30 MPa for steel over several mm^ in area. 
In the blind-hole drilling method a special rosette strain gage is 
bonded to the test specimen. A small hole (typically 1.5 to 3 mm in 
both diameter and depth) is then drilled through the center of the gage 
and the resulting change in strain around the hole due to relaxation is 
measured. Residual stresses are obtained from the strains by analytic 
methods. Blind-hole drilling is a semi-destructive technique and it 
measures the average strain within the top surface layers at a very 
specific location. This technique is limited to thick samples and 
cannot be used for thin silicon plates since the drilling procedure can 
noticeably impair the structural integrity of the part being tested. 
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Both techniques are unsatisfactory for measuring residual stresses 
because only stresses within the top few millimeters at specific 
locations are measured and full field, through the tickness residual 
stresses are usually desired. Moreover, considering the additional 
limitations imposed by thin plates it can be concluded that these 
requirements can present formidable obstacles in determining residual 
stresses . 

In this paper, a laser interferometry technique [14,15] for measuring 
deflections over large spatial dimensions is described. These 
measurements, when combined with an analysis for predicting deflections 
in linear-elastic plates, lead to the residual stresses. This paper 
outlines the procedure for obtaining the residual stresses from the 
experimental measurements in thin circular plates. 

Description of the Technique 

The technique largely depends on the ability to measure deflections of 
thin plates when the applied loads are known. Figure 1 shows a 
schematic of the supporting and loading geometry. A circular thin plate 
is simply supported along the edge end loaded at the center. The true 
deflections can be measured by laser interferometry; the experimental 
arrangement of which is shown in Fig. 2. The coherent light produced by 
a helium-neon laser is diverged and collimated by an optical system. 
The uniformly expanded beam (10 cm in diameter) emerging from the 
parabolic mirror is reflected onto the specimen by a broadband 
beamsplitter inclined at 45° to the optical axis. The return beam from 
the polished surface of the specimen is converged by a convex lens after 
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passing through the optical flat and the beamsplitter. Finally, the 
converged light beam carrying the optical information (interference 
pattern) is recorded by a camera. 


The constructive and destructive interference of the monochromatic light 
source occurs in the following way. The monochromatic light is directed 
perpendicular to the specimen and reflected from the bottom surface of 
the optical flat and the top surface of the specimen and viewed in 
reflection. Destructive interference (dark fringes) occurs when the 
optical path difference (2P) is an odd number of half wavelengths V2 



n = 1,2,3,... 


Constructive 


interference (light fringes) occurs when 


( 1 ) 


2 p = n A 


n 1 ,2 , 3 , . • . 


The interference pattern can be analyzed to yield full-field plate 
deflection to an accuracy of half wavelength. 


Determination of Res idual Stre s s e s from De flections 


Consider a plate element such as shown in Fig. 3. Taking the sum of the 
projections in the radial direction of the in-plane forces acting on the 
element, we obtain 

dN r ' (2) 

N r - N t + JT = ° 

where N f = c^h and = cr t h are the radial and tangential components 
of the stress residuals, 0 r and 0 t 


the average residual stresses and h 
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the plate thickness. Also, as given in reference [16] the governing 
plate equation can be expressed in the following form 





1 du)\ 

7 7J 



+ 


p 

2TTr 


(3) 


where 


Eh“ 


D = 


12(1 - vi ) 


E = Young's modulus 
Vi = Poisson's ratio 
w = deflection 

r = radial distance measured from point of loading 
Solving for we get 



(4) 


Hence the residual stresses may be obtained if the deflection is known 
as a function of r. 


The experimentally determined deflection (l, ex p obtained from the 
interference pattern may be represented in the simplest functional form 
as a third degree polynomial 

(a = Ar 3 + Br 2 + Cr + F 

in the vicinity of any point Q (see Fig. 4). The constants A, B, C 
and F may be found given at ( i = 1,2, 3,4) as 
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Once A, B, C and F are found the derivatives of the experimenal 
deflection curve at point Q can be evaluated as 


do) -) 

dr Q = 3Ar Q + 2Br Q + C 


d 2 cx) 


= 6 Ar n + 2 B 
Q 0 




6A 


( 6 ) 


Equations (6) can be substituted into Eq. (4) to give the experimental 
N r at point Q as 


D 9A- — - 


2ti rO 


N r " 


3 Ar q + 2Btq + C 


( 7 ) 


This procedure can be repeated (n-3) times if there are n experimental 
data points. The values of can then be plotted versus r. After 
plotting the experimental values of along a radial line, Eq . (2) 
can be used in conjunction with this plot to give the experimental 


distribution. 
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Experiments 

The technique described above was used to determine residual stresses in 
a single crystal, 0.575 mm thick and 6.25 cm in diameter Czochralski 
silicon wafer. The optical system presented in Fig. 2 and a fixture 
capable of testing circular wafers under axisymmetric loading were 
built. The specimen was simply supported with a supporting edge 
diameter of 89 mm (2a). The interference pattern before loading and 
under the applied load was recorded with a 4 x 5 camera. A typical 
interference pattern for the silicon wafer is shown in Fig. 5. Each 
fringe spacing represents an incremental deflection of 315 nm which is 
half the wavelength of the lle-Ne laser. A nikon comparator ( 6C-2 ) was 
used to determine the fringe locations with high accuracy. Values of 
fringe order and corresponding centroidal distances were recorded to 
generate surface profiles along eight radial lines. Local deflections 
were found and fed into a computer program to evaluate stress residuals 
along these lines according to the analysis described in the previous 
section. 

Results and Discuss ion 

A data set including fringe orders, radial distances, deflections and 
residual stress values are presented in Table I. The radial distances 
were measured from the point of load application. The central 
deflection of the plate, w ma ^ , was measured independently. The fringe 
orders were numbered such that the first contour around the point of 
loading was of zero order. The deflections were then found as: 
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X 

w A = w_ - n A — 

Q max Q 2 

where 

Wq = deflection at point Q 
n Q = fringe order at point Q 
^ = wavelength (628 nm) 

The combined resolution of the optical setup and the comparator in 
determining the fringe locations was 1 pm. 

A topographic scheme was used to plot the constant value contour line 
systems for the average residual stress which is shown in Fig. 6. 

Only the central region of the circular wafer was analyzed which 
represents an area of 28 cm . The sensitivity of this integrated 
technique in determining residual stresses is found to vary between 
lOOKPa and 500KPa. Applied flexural load, plate curvature, fringe 
density and optical alignment can be listed among the parameters 
effecting the uncertainty. A detailed sensitivity analysis for 
different plate geometries will be the subject of a subsequent paper. 
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Fig . 1 

Simply supported circular wafer 
under centric loading 


Fig. 2 


Schematic of the experimental setup 



Fig . 3 Fig . 4 

Definition of the stress Presentation of data points around 

residuals on a plate element point Q where residual stresses 

in polar coordinates are evaluated 
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Fig. 5 Fig. 6 

Typical interference pattern Distribution of the average 

residual stress 0_ 


TABLE I 

Typical Data Set 


1 1 

1 Fringe 1 

1 

Radial ! 

Deflection 

1 

I Avg . Residual Stress 

1 Order 1 

Distance 1 

w ( m) 

1 a r 


5 

0 .009674 

0.000051322 


6 

0 .009983 

0.000051008 

5 . 5 MPa 1 

7 

0.010287 

0.000050694 


8 

0.010587 

0.000050380 

(at Tq = 0.01m) 1 

37 

0.014641 

0.000045670 


38 

0.014890 

0.000045356 

-2.66 MPa 1 

39 

0.015136 

0.000045042 


40 

0.015377 

0.000044728 

(at Tq = 0 .015m) I 








